Methods of obtaining ophthalmic lenses providing the eye with reduced aberrations 
Field of invention 

The present invention relates to methods of designing ophthalmic lenses that 
provide the eye with reduced aberrations, as well as lenses capable of providing such 
visual improvements. 

Background of the invention 

It is presently discussed that the visual quality of eyes having an implanted 
intraocular lens (IOL) is comparable with normal eyes in a population of the same age. 
Consequently, a 70 year old cataract patient can only expect to obtain the visual quality 
of a non-cataracteous person of the same age after surgical implantation of an intraocular 
lens, although such lenses objectively have been regarded as optically superior to the 
natural crystalline lens. This result is likely to be explained by the fact that present IOLs 
are not adapted to compensate for defects of the optical system of the human eye, namely 
optical aberrations. Age-related defects of the eye have recently been investigated and it 
is found that contrast sensitivity significantly declines in subjects older than 50 years. 
These results seem to comply with the above-mentioned discussion, since the contrast 
sensitivity measurements indicate that individuals having undergone cataract surgery with 
lens implantation lens will not obtain a better contrast sensitivity than non-cataracteous 
persons of an average age of about 60 to 70 years. 

Even if intraocular lenses aimed at substituting the defective cataract lens and 
other ophthalmic lenses, such as conventional contact lenses, have been developed with 
excellent optical quality as isolated elements, it is obvious that they fail to correct for a 
number of aberration phenomena of the eye including age-related aberration defects. 

US Patent No. 5,777,719 (Williams et al.) discloses a method and an apparatus for 
accurately measuring higher order aberrations of the eye as an optical system with 
wavefront analysis. By using a Hartmann-Shack wavefront sensor, it is possible to 
measure higher order aberrations of the eye and use such data to find compensation for 
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these aberrations and thereby obtain sufficient information for the design of an optical 
lens, which can provide a highly improved optical performance. The Hartmann-Shack 
sensor provides means for analyzing light reflected from a point on the retina of the eye 
of a subject. The wavefront in the plane of the pupil is recreated in the plane of the lenslet 
5 array of the Hartmann-Shack sensor. Each lenslet in the array is used to form an aerial 
image of the retinal point source on a CCD camera located at the focal plane of the array. 
The wave aberration of the eye, in the form resulting from a point source produced on 
the retina by a laser beam, displaces each spot by an amount proportional to the local 
slope of the wavefront at each of the lenslets. The output from the CCD camera is sent to 
10 a computer, which then performs calculations to fit slope data to the first derivatives of 
66 Zernike polynomials. From these calculations, coefficients for weighting the Zernike 
polynomials are obtained. The sum of the weighted Zernike polynomials represents a 
reconstructed wavefront distorted by the aberrations of the eye as an optical system. The 
individual Zernike polynomial terms will then represent different modes of aberration. 
15 US Patent No. 5,050,981 (Roffman) discloses another method for designing a lens 

by calculating modulation transfer functions from tracing a large number of rays through 
the lens-eye system and evaluating the distribution density of the rays in the image 
position. This is repeatedly performed by varying at least one lens surface until a lens is 
found which results in a sharp focus and a maximum modulation transfer function. 
20 US Patent No. 6,224,21 1 (Gordon) describes a method of improving the visual 

acuity of the human eye by successively fitting aspheric lenses to the cornea and thereby 
finding a lens that can reduce spherical aberration of the whole individual eye. 

The methods referred to above for designing are suitable for the design of contact 
lenses or other correction lenses for the phakic eye which can be perfected to compensate 
25 for the aberration of the whole eye system. However, to provide improved intraocular 
lenses aimed to replace the natural crystalline lens, it would be necessary to consider the 
aberrations of the individual parts of the eye. 

US Patent No. 6,050,687 (Bille et al) refers to a method wherein the refractive 
properties of the eye are measured and wherein consideration is taken to the contribution 
30 of the individual surfaces of the eye to the total wavefront aberrations. The method 



3 



described herein particularly aims at analyzing the topography of the posterior corneal 
surface in order to improve refractive correction techniques. 

There has recently been a focus on studying the aberrations of the eye, including a 
number of studies of the development of these aberrations as a function of age. In two 
5 particular studies, the development of the components of the eye were examined 
separately, leading to the conclusion that the optical aberrations of the individual 
components of younger eyes cancel each other out, see Optical Letters, 1998, Vol. 
23(21), pp.1713-1715 and IOVS, 2000, Vol.41(4), 545. The article of S. Patel et al in 
Refractive & Corneal Surgery, 1993, Vol. 9, pages 173-181 discloses the asphericity of 
10 posterior corneal surfaces. It is suggested that the corneal data can be used together with 
other ocular parameters to predict the power and the asphericity of an intraocular lens 
with the purpose of maximizing the optical performances of the future pseudophakic eye. 
Furthermore, it was also recently observed by Antonio Guirao and Pablo Artal in IOVS, 
1999, Vol. 40(4), S535 that the shape of the cornea changes with age and becomes more 
15 spherical. These studies indicate that the cornea in the subjects provides a positive 
spherical aberration, which increases slightly with the age. On the other hand, the 
rotationally symmetric aberration of the anterior corneal surface does not seem to be 
different between younger and older eye according to results found by T Oshika et al in 
Investigative Ophthalmology and Visual Science, 1999, Vol. 40, pp. 135 1-1 355. In Vision 
20 Research, 1998, 38(2), pp. 209-229, A Glasser et al. investigated the spherical aberration 
of natural crystalline lenses from eyes obtained from an eye bank after the cornea has 
been removed. According to the laser scanner optical method used herein it was found 
that the spherical aberration from an older lens (66 years) shows positive spherical 
aberration, whereas a 1 0-year-old lens shows negative spherical aberration. In addition, 
25 Vision Research, 2001 , 41 , pp.235-243 (G Smith et al) discloses that the natural 

crystalline lens appears to have negative spherical aberration when in the relaxed state. 
Smith et al suggest that because older eyes have a larger aberration, it is likely that the 
spherical aberration of the crystalline lens becomes less negative with age. 

In Ophthal. Physiol. Opt, 1991, Vol. 11, pp. 137-143 (DA Atchison) it is 
30 discussed how to reduce spherical aberrations in intraocular lenses by aspherizing the 
lens surface. The methods outlined by Atchison are based on geometric transferring 



calculations, which do not consider diffraction effects and any variations in refractive 
index along the ray path in inhomogeneous elements. These calculations will lead to 
errors close to the diffraction limit. Also in WO 98/31299 (Technomed) a ray tracing 
method is outlined according to which the refraction of the cornea is attempted to be 
considered for the design of an intraocular lens.In view of the foregoing, it is apparent 
that there is a need for ophthalmic lenses that are better adapted or compensated to the 
aberrations of the individual surfaces of the eye and are capable of better correcting 
aberrations other than defocus and astigmatism, as provided by conventional ophthalmic 
lenses. 

Description of the invention 

It is an object of the invention to provide for methods that result in obtaining an 
ophthalmic lens, which provides the eye with reduced aberrations. 

It is another object of the invention to provide methods of obtaining an intraocular 
lens capable of reducing the aberration of the eye after its implantation into the eye. 

It is a further object to provide for methods of obtaining an intraocular lens 
capable of compensating the aberrations resulting from optical irregularities in the 
corneal surfaces. 

It is a still further object of the present invention to provide an intraocular lens 
which is capable of restoring a wavefront deviating from sphericity into a substantially 
more spherical wavefront. 

It is a also an object of the present invention to provide an intraocular lens which 
is capable of correcting for mean optical irregularities and imperfections found in a 
particular group of people and thereby provide a lens with improved optical performance 
for an individual belonging to the same group. 

The present invention generally relates to an ophthalmic lens and to methods of 
obtaining said ophthalmic lens that is capable of reducing the aberrations of the eye. By 
aberrations in this context is meant wavefront aberrations. This is based on the 
understanding that a converging wavefront must be perfectly spherical to form a point 
image, i.e. if a perfect image shall be formed on the retina of the eye, the wavefront 
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having passed the optical surfaces of the eye, such as the cornea and a natural or artificial 
lens, must be perfectly spherical. An aberrated image will be formed if the wavefront 
deviates from being spherical. . In this context the term nonspherical surface will refer to 
rotationally symmetric, asymmetric and/or irregular surfaces, i.e. all surfaces differing 
5 from a sphere. The wavefront aberrations can be expressed in mathematical terms in 
accordance with different approximate models as is explained in textbook references, 
such as M.R. Freeman, Optics, Tenth Edition, 1990. 

In a first embodiment, the present invention is directed to a method of designing 
an intraocular lens capable of reducing aberrations of the eye after its implantation. The 

10 method comprises a first step of characterizing at least one corneal surface as a 

mathematical model and by employing the mathematical model calculating the resulting 
aberrations of the corneal surface. An expression of the corneal aberrations is thereby 
obtained, i.e. the wavefront aberrations of a spherical wavefront having passed such a 
corneal surface. Dependent on the selected mathematical model different routes to 

15 calculate the corneal aberrations can be taken. Preferably, the corneal surface is 

characterized as a mathematical model in terms of a conoid of rotation or in terms of 
polynomials or a combination thereof. More preferably, the corneal surface is 
characterized in terms of a linear combination of polynomials. The second step of the 
method is to select the power of the intraocular lens, which is done according to 

20 conventional methods for the specific need of optical correction of the eye, for example 
the method described in US Patent No. 5,968,095 From the information of steps one and 
two an intraocular lens is modeled, such that a wavefront from an optical system 
comprising said lens and corneal model obtains reduced aberrations. The optical system 
considered when modeling the lens typically includes the cornea and said lens, but in the 

25 specific case it can also include other optical elements including the lenses of spectacles, 
or an artificial correction lens, such as a contact lens, a corneal inlay implant or an 
implantable correction lens depending on the individual situation. 

Modeling the lens involves selection of one or several lens parameters in a system 
which contributes to determine the lens shape of a given, pre-selected refractive power. 

30 This typically involves the selection of the anterior radius and surface shape^ posterior 
radius and surface shape, the lens thickness and the refractive index of the lens. In 
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practical terms, the lens modeling can be performed with data based on a conventional 
spherical lens, such as the CeeOn® lenses from Pharmacia Corp., as exemplified with the 
CeeOn® Edge (Model 91 1). In such a case, it is preferred to deviate as little as possible 
from an already clinically approved model. For this reason, it may be preferred to 
5 maintain pre-determined values of the central radii of the lens, its thickness and refractive 
index, while selecting a different shape of the anterior and/or posterior surface, thus 
providing one or both of these surfaces to have an nonspherical shape. According to an 
alternative of the inventive method, the spherical anterior surface of the conventional 
starting lens is modeled by selecting a suitable aspheric component. Preferably the lens 
10 has at least one surface described as a nonsphere or other conoid of rotation. Designing 
nonspherical surfaces of lenses is a well-known technique and can be performed 
according to different principles and the description of such surfaces is explained in more 
detail in our parallel Swedish Patent Application 000061 1-4 to which is given reference. 
The inventive method can be further developed by comparing wavefront 
1 5 aberrations of an optical system comprising the lens and the model of the average cornea 
with the wavefront aberrations of the average cornea and evaluating if a sufficient 
reduction in wavefront aberrations is obtained. Suitable variable parameters are found 
among the above-mentioned physical parameters of the lens, which can be altered so as to 
find a lens model, which deviates sufficiently from being a spherical lens to compensate 
20 for the corneal aberrations. 

The characterization of at least one corneal surface as a mathematical model and 
thereby establishing a corneal model expressing the corneal wavefront aberrations is 
preferably performed by direct corneal surface measurements according to well-known 
topographical measurement methods which serve to express the surface irregularities of 
25 the cornea in a quantifiable model that can be used with the inventive method. Corneal 
measurements for this purpose can be performed by the ORBSCAN® videokeratograph, 
as available from Orbtech, , or by corneal topography methods, such as EyeSys® from 
Premier Laser Systems. Preferably, at least the front corneal surface is measured and 
more preferably both front and rear corneal surfaces are measured and characterized and 
30 expressed together in resulting wavefront aberration terms, such as a linear combination 
of polynomials which represent the total corneal wavefront aberrations. According to one 
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important aspect of the present invention, characterization of corneas is conducted on a 
selected population with the purpose of expressing an average of corneal wavefront 
aberrations and designing a lens from such averaged aberrations. Average corneal 
wavefront aberration terms of the population can then be calculated, for example as an 
5 average linear combination of polynomials and used in the lens design method. This 
aspect includes selecting different relevant populations, for example in age groups, to 
generate suitable average corneal surfaces. Advantageously, lenses can thereby be 
provided which are adapted to an average cornea of a population relevant for an 
individual elected to undergo cataract surgery or refractive correction surgery including 
10 implantation of an IOL or corneal inlays. The patient will thereby obtain a lens that gives 
the eye substantially less aberrations when compared to a conventional spherical lens. 

Preferably, the mentioned corneal measurements also include the measurement of 
the corneal refractive power. The power of the cornea and the axial eye length are 
typically considered for the selection of the lens power in the inventive design method. 
15 Also preferably, the wavefront aberrations herein are expressed as a linear 

combination of polynomials and the optical system comprising the corneal model and 
modeled intraocular lens provides for a wavefront having obtained a substantial reduction 
in aberrations, as expressed by one or more such polynomial terms. In the art of optics, 
several types of polynomials are available to skilled persons for describing aberrations. 
20 Suitably, the polynomials are Seidel or Zernike polynomials. According to the present 
invention Zernike polynomials preferably are employed. 

The technique of employing Zernike terms to describe wavefront aberrations 
originating from optical surfaces deviating from being perfectly spherical is a state of the 
art technique and can be employed for example with a Hartmann-Shack sensor as 
25 outlined in J. Opt. Soc. Am., 1994, Vol. 1 1(7), pp. 1949-57. It is also well established 
among optical practitioners that the different Zernike terms signify different aberration 
phenomena including defocus, astigmatism, coma and spherical aberration up to higher 
aberrations. In an embodiment of the present method, the corneal surface measurement 
results in that a corneal surface is expressed as a linear combination of the first 1 5 
30 Zernike polynomials. By means of a raytracing method, the Zernike description can be 
transformed to a resulting wavefront (as described in Equation (1)), wherein Zj is the i-th 
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Zernike term and aj is the weighting coefficient for this term. Zernike polynomials are a 
set of complete orthogonal polynomials defined on a unit circle. Below, Table 1 shows 
the first 15 Zernike terms and the aberrations each term signifies. 

5 = jf>,Z, (1) 

In equation (l),p and 0 represent the normalized radius and the azimuth angle, 
respectively. 

10 Table 1 
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Conventional optical correction with intraocular lenses will only comply with the 
fourth term of an optical system comprising the eye with an implanted lens. Glasses, 
contact lenses and some special intra ocular lenses provided with correction for 
astigmatism can further comply with terms five and six and substantially reducing 
Zemike polynomials referring to astigmatism. 

The inventive method further includes to calculate the wavefront aberrations 
resulting from an optical system comprising said modeled intraocular lens and cornea and 
expressing it in a linear combination of polynomials and to determine if the intraocular 
lens has provided sufficient reduction in wavefront aberrations. If the reduction in 
wavefront aberrations is found to be insufficient, the lens will be re-modeled until one or 
several of the polynomial terms are sufficiently reduced. Remodeling the lens means that 
at least one lens design parameter is changed. These include the anterior surface shape 
and central radius, the posterior surface shape and central radius, the thickness of the lens 
and its refractive index. Typically, such remodeling includes changing the shape of a lens 
surface so it deviates from being a spherical. There are several tools available in lens 
design that are useful to employ with the design method, such as OSLO version 5 see 
Program Reference, Chapter 4, Sinclair Optics 1996. The format of the Zernike 
polynomials associated with this application are listed in Table 1 . 

According to a preferred aspect of the first embodiment, the inventive method 
comprises expressing at least one corneal surface as a linear combination of Zemike 
polynomials and thereby determining the resulting corneal wavefront Zemike 
coefficients, i.e. the coefficient of each of the individual Zemike polynomials that is 
selected for consideration. The lens is then modeled so that an optical system comprising 
of said model lens and cornea provides a wavefront having a sufficient reduction of 
selected Zernike coefficients. The method can optionally be refined with the further steps 
of calculating the Zemike coefficients of the Zemike polynomials representing a 
wavefront resulting from an optical system comprising the modeled intraocular lens and 
cornea and determining if the lens has provided a sufficient reduction of the wavefront 
Zemike coefficients of the optical system of cornea and lens; and optionally re-modeling 
said lens until a sufficient reduction in said coefficients is obtained. Preferably, in this 
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aspect the method considers Zernike polynomials up to the 4th order and aims to 
sufficiently reduce Zernike coefficients referring to spherical aberration and/or 
astigmatism terms. It is particularly preferable to sufficiently reduce the 1 1th Zernike 
coefficient of a wavefront front from an optical system comprising cornea and said 
modeled intraocular lens, so as to obtain an eye sufficiently free from spherical 
aberration. Alternatively, the design method can also include reducing higher order 
aberrations and thereby aiming to reduce Zernike coefficients of higher order aberration 
terms than the 4 th order. 

When designing lenses based on corneal characteristics from a selected 
population, preferably the corneal surfaces of each individual are expressed in Zernike 
polynomials describing the surface topography and therefrom the Zernike coefficients of 
the wavefront aberration are determined. From these results average Zernike wavefront 
aberration coefficients are calculated and employed in the design method, aiming at a 
sufficient reduction of selected such coefficients. In an alternative method according to 
the invention, average values of the Zernike polynomials describing the surface 
topography are instead calculated and employed in the design method. It is to be 
understood that the resulting lenses arriving from a design method based on average 
values from a large population have the purpose of substantially improving visual quality 
for all users. A lens having a total elimination of a wavefront aberration term based on an 
average value may consequently be less desirable and leave certain individuals with an 
inferior vision than with a conventional lens. For this reason, it can be suitable to reduce 
the selected Zernike coefficients only to certain degree of the average value. 

According to another approach of the inventive design method, corneal 
characteristics of a selected population and the resulting linear combination of 
polynomials, e.g. Zernike polynomials, expressing each individual corneal aberration can 
be compared in terms of coefficient values. From this result, a suitable value of the 
coefficients is selected and employed in the inventive design method for a suitable lens. 
In a selected population having aberrations of the same sign such a coefficient value can 
typically be the lowest value within the selected population and the lens designed from 
this value would thereby provide improved visual quality for all individuals in the group 
compared to a conventional lens. One embodiment of the method comprises selecting a 
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representative group of patients and collecting corneal topographic data for each subject 
in the group. The method comprises further transferring said data to terms representing 
the corneal surface shape of each subject for a preset aperture size representing the pupil 
diameter. Thereafter a mean value of at least one corneal surface shape term of said group 
is calculated, so as to obtain at least one mean corneal surface shape term. Alternatively 
or complementary a mean value of at least one to the cornea corresponding corneal 
wavefront aberration term can be calculated. The corneal wavefront aberration terms are 
obtained by transforming corresponding corneal surface shape terms using a raytrace 
procedure. From said at least one mean corneal surface shape term or from said at least 
one mean corneal wavefront aberration term an ophthalmic lens capable of reducing said 
at least one mean wavefront aberration term of the optical system comprising cornea and 
lens is designed. 

In one preferred embodiment of the invention the method further comprises 
designing an average corneal model for the group of people from the calculated at least 
one mean corneal surface shape term or from the at least one mean corneal wavefront 
aberration term. It also comprises checking that the designed ophthalmic lens 
compensates correctly for the at least one mean aberration term. This is done by 
measuring these specific aberration terms of a wavefront having traveled through the 
model average cornea and the lens. The lens is redesigned if said at least one aberration 
term has not been sufficiently reduced in the measured wavefront. 

Preferably one or more surface descriptive (asphericity describing) constants are 
calculated for the lens to be designed from the mean corneal surface shape term or from 
the mean corneal wavefront aberration terms for a predetermined radius. The spherical 
radius is determined by the refractive power of the lens. 

The corneal surfaces are preferably characterized as mathematical models and the 
resulting aberrations of the corneal surfaces are calculated by employing the 
mathematical models and raytracing techniques. An expression of the corneal wavefront 
aberrations is thereby obtained, i.e. the wavefront aberrations of a wavefront having 
passed such a corneal surface. Dependent on the selected mathematical model different 
routes to calculate the corneal wavefront aberrations can be taken. Preferably, the corneal 
surfaces are characterized as mathematical models in terms of a conoid of rotation or in 
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terms of polynomials or a combination thereof. More preferably, the corneal surfaces are 
characterized in terms of linear combinations of polynomials. 

In one embodiment of the invention, the at least one nonspheric surface of the lens 
is designed such that the lens, in the context of the eye, provides to a passing wavefront at 
least one wavefront aberration term having substantially the same value but with opposite 
sign to a mean value of the same aberration term obtained from corneal measurements of 
a selected group of people, to which said patient is categorized. Hereby a wavefront 
arriving from the cornea of the patient's eye obtains a reduction in said at least one 
aberration term provided by the cornea after passing said lens. The used expression c in 
the context of the eye' can mean both in the real eye and in a model of an eye. 
In a specific embodiment of the invention, the wavefront obtains reduced aberration 
terms expressed in rotationally symmetric Zemike terms up to the fourth order. For this 
purpose, the surface of the ophthalmic lens is designed to reduce a positive spherical 
aberration term of a passing wavefront. The consequence of this is that if the cornea is a 
perfect lens and thus not will give rise to any wavefront aberration terms the ophthalmic 
lens will provide the optical system comprising the cornea and the ophthalmic lens with a 
negative wavefront spherical aberration term. In this text positive spherical aberration is 
defined such that a spherical surface with positive power produces positive spherical 
aberration. Preferably the lens is adapted to compensate for spherical aberration and more 
preferably it is adapted to compensate for at least one term of a Zernike polynomial 
representing the aberration of a wavefront, preferably at least the 1 1 th Zemike term, see 
Table 1. 

The selected groups of people could for example be a group of people belonging 
to a specific age interval, a group of people who will undergo a cataract surgical 
operation or a group of people who have undergone corneal surgery including but not 
limited to LASIK (laser in situ keratomileusis), RK (radial keratotomy) or PRK 
(photorefractive keratotomy). The group could also be a group of people who have a 
specific ocular disease or people who have a specific ocular optical defect. 

The lens is also suitably provided with an optical power. This is done according 
to conventional methods for the specific need of optical correction of the eye. Preferably 
the refractive power of the lens is less than or equal to 30 diopters. An optical system 
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considered when modeling the lens to compensate for aberrations typically includes the 
average cornea and said lens, but in the specific case it can also include other optical 
elements including the lenses of spectacles, or an artificial correction lens, such as a 
contact lens, a corneal inlay or an implantable correction lens depending on the individual 
situation. 

In an especially preferred embodiment the ophthalmic lens is designed for people 
who will undergo a cataract surgery. In this case it is has been shown that the average 
cornea from such a population is represented by a prolate surface following the formula: 



wherein, 

(i) the conical constant cc has a value ranging between -1 and 0 

(ii) R is the central lens radius and 

(iii) ad and ae are aspheric polynomial coefficients additional to the conical constant. 
In these studies the conic constant of the prolate surface ranges between about — 

0.05 for an aperture size (pupillary diameter) of 4 mm to about -0. 1 8 for an aperture size 
of 7 mm. Accordingly an ophthalmic lens suitable to improve visual quality by reducing 
at least spherical aberration for a cataract patient based on an average corneal value will 
have a prolate surface following the formula above. Since cornea generally produces a 
positive spherical aberration to a wavefront in the eye, an ophthalmic lens for 
implantation into the eye will have negative spherical aberration terms while following 
the mentioned prolate curve. As will discussed in more detail in the exemplifying part of 
the specification, it has been found that an intraocular lens that can correct for 100% of a 
mean spherical aberration has a conical constant (cc) with a value of less than 0 
(representing a modified conoid surface), with an exact value dependent on the design 
pupillary diameter and the selected refractive power. For example, a 6 mm diameter 
aperture will provide a 22 diopter lens with conical constant value of about — 1.03. 
In this embodiment, the ophthalmic lens is designed to balance the spherical aberration of 
a cornea that has a Zemike polynomial coefficient representing spherical aberration of the 
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wavefront aberration with a value in the interval from 0.000156 mm to 0.001948 mm for 
a 3mm aperture radius, 0.000036 mm to 0.000448 mm for a 2 mm aperture radius, 
0.0001039 mm to 0.0009359 mm for a 2,5 mm aperture radius and 0.000194 mm to 
0.00365 mm for a 3,5 mm aperture radius using polynomials expressed in OSLO format. 
These values were calculated for a model cornea having a single surface with a refractive 
index of the cornea of 1 .3375. It is possible to use optically equivalent model formats of 
the cornea without departing from the scope of the invention. For example multiple 
surface corneas or corneas with different refractive indices could be used. The lower 
values in the intervals are here equal to the measured average value for that specific 
aperture radius minus one standard deviation. The higher values are equal to the 
measured average value for each specific aperture radius plus three standard deviations. 
The used average values and standard deviations are shown in tables 8,9,10 and 11. The 
reason for selecting only minus one SD ( = Standard Deviation) while selecting plus three 
SD is that in this embodiment it is convenient to only compensate for positive corneal 
spherical aberration and more than minus one SD added to the average value would give 
a negative corneal spherical aberration. 

According to one embodiment of the invention the method further comprises the 
steps of measuring the at least one wavefront aberration term of one specific patient's 
cornea and determining if the selected group corresponding to this patient is 
representative for this specific patient. If this is the case the selected lens is implanted and 
if this is not the case a lens from another group is implanted or an individual lens for this 
patient is designed using this patients corneal description as a design cornea. These 
method steps are preferred since then patients with extreme aberration values of their 
cornea can be given special treatments. 

According to another embodiment, the present invention is directed to the 
selection of an intraocular lens of refractive power, suitable for the desired optical 
correction that the patient needs, from a plurality of lenses having the same power but 
different aberrations. The selection method is similarly conducted to what has been 
described with the design method and involves the characterizing of at least one corneal 
surface with a mathematical model by means of which the aberrations of the corneal 
surface is calculated. The optical system of the selected lens and the corneal model is 
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then evaluated so as to consider if sufficient reduction in aberrations is accomplished by 
calculating the aberrations of a wavefront arriving from such a system. If an insufficient 
correction is found a new lens is selected, having the same power, but different 
aberrations. The mathematical models employed herein are similar to those described 
above and the same characterization methods of the corneal surfaces can be employed. 

Preferably, the aberrations determined in the selection are expressed as linear 
combinations of Zernike polynomials and the Zernike coefficients of the resulting optical 
system comprising the model cornea and the selected lens are calculated. From the 
coefficient values of the system, it can be determined if the intraocular lens has 
sufficiently balanced the corneal aberration terms, as described by the Zernike 
coefficients of the optical system. If no sufficient reduction of the desired individual 
coefficients are found these steps can be iteratively repeated by selecting a new lens of 
the same power but with different aberrations, until a lens capable of sufficiently 
reducing the aberrations of the optical system is found. Preferably at least 1 5 Zernike 
polynomials up to the 4 th order are determined. If it is regarded as sufficient to correct for 
spherical aberration, only the spherical aberration terms of the Zernike polynomials for 
the optical system of cornea and intraocular lens are corrected. It is to be understood that 
the intraocular lens shall be selected so a selection of these terms become sufficiently 
small for the optical system comprising lens and cornea. In accordance with the present 
invention, the 1 1 th Zernike coefficient, an, can be substantially eliminated or brought 
sufficiently close to zero. This is a prerequisite to obtain an intraocular lens that 
sufficiently reduces the spherical aberration of the eye. The inventive method can be 
employed to correct for other types of aberrations than spherical aberration by 
considering other Zernike coefficients in an identical manner, for example those 
signifying astigmatism, coma and higher order aberrations. Also higher order aberrations 
can be corrected dependent on the number of Zernike polynomials elected to be a part of 
the modeling, in which case a lens can be selected capable of correcting for higher order 
aberrations than the 4 th order. 

According to one important aspect, the selection method involves selecting lenses 
from a kit of lenses having lenses with a range of power and a plurality of lenses within 
each power having different aberrations. In one example the lenses within each power 
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have anterior surfaces with different aspherical components. If a first lens does not 
exhibit sufficient reduction in aberration, as expressed in suitable Zernike coefficients, 
then a new lens of the same power, but with a different surface (aspheric component) is 
selected. The selection method can if necessary be iteratively repeated until the best lens 
5 is found or the studied aberration terms are reduced below a significant borderline value. 
In practice, the Zernike terms obtained from the corneal examination will be directly 
obtained by the ophthalmic surgeon and by means of an algorithm they will be compared 
to known Zernike terms of the lenses in the kit. From this comparison the most suitable 
lens in the kit can be found and implanted. Alternatively, the method can be conducted 
10 before cataract surgery and data from the corneal estimation is sent to a lens manufacturer 
for production of an individually tailored lens. 

The present invention further pertains to an intraocular lens having at least one 
nonspherical surface capable of transferring a wavefront having passed through the 
cornea of the eye into a substantially spherical wavefront with its center at the retina of 
15 the eye. Preferably, the wavefront is substantially spherical with respect to aberration 
terms expressed in rotationally symmetric Zernike terms up to the fourth order. 

In accordance with an especially preferred embodiment, the invention relates to 
an intraocular lens, which has at least one surface, when expressed as a linear 
combination of Zernike polynomial terms using the normalized format, that has a 
20 negative 1 1 th term of the fourth order with a Zernike coefficient ai i that that can balance a 
positive such term of the cornea n to obtain sufficient reduction of the spherical 
aberration of the eye after implantation. In one aspect of this embodiment, the Zernike 
coefficient an of the lens is determined so as to compensate for an average value resulting 
from a sufficient number of estimations of the Zernike coefficient an in several corneas. 
25 In another aspect, the Zernike coefficient ai \ is determined to compensate for the 

individual corneal coefficient of one patient. The lens can accordingly be tailored for an 
individual with high precision. 

The invention further relates to another method of providing a patient with an 
intraocular lens, which at least partly compensates for the aberrations of the eye. This 
30 method comprises removing the natural lens from the eye. Surgically removing the 
impaired lens can be performed by using a conventional phacoemulsification method. 
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The method further comprises measuring the aberrations of the aphakic eye, not 
comprising a lens, by using a wavefront sensor. Suitable methods for wavefront 
measurements are found in J.Opt.Soc.Am., 1994, Vol. 1 1(7), pp. 1949-57 by Liang et. al. 
Furthermore, the method comprises selecting from a kit of lenses a lens that at least partly 
compensates for the measured aberrations and implanting said lens into the eye. The kit 
of lenses comprises lenses of different power and different aberrations and finding the 
most suitable lens can be performed in a manner as earlier discussed. Alternatively, an 
individually designed lens for the patient can be designed based on the wavefront analysis 
of the aphakic eye for subsequent implantation. This method is advantageous, since no 
topographical measurements of the cornea are and the whole cornea, including the front 
and back surfaces, is automatically considered. 

The lenses according to the present invention can be manufactured with 
conventional methods. In one embodiment they are made from soft, resilient material, 
such as silicones or hydrogels. Examples of such materials suitable for foldable 
intraocular lenses are found in US patent No. 5,444,106 or in US Patent No. 5,236,970. 
Manufacturing of nonspherical silicone lenses or other foldable lenses can be performed 
according to US Patent No. 6,007,747. Alternatively, the lenses according to the present 
invention can be made of a more rigid material, such as polymethylmethacrylate. The 
skilled person can readily identify alternative materials and manufacturing methods, 
which will be suitable to employ to produce the inventive aberration reducing lenses. 
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Detailed description of the invention 

Fig.l shows a comparison of the an ("Zl 1") Zernike coefficient values for 10 subjects if 
implanted with CeeOn® 911 lenses and the inventive averaged ("Zl 1") lens. 
Fig. 2 shows modeled visual acuities of the test subjects with CeeOn® 911 lenses and the 
inventive averaged ("Zl 1") lenses. 

Fig. 3 and Fig. 4 show modulation transfer function comparisons between CeeOn® 91 1 
lenses and the inventive averaged ("Zl 1") lenses 

Fig. 5 shows visual acuity plotted as a function of the astigmatism of the lenses according 
to the model lenses according to the invention. 

Fig. 6 shows the best corrected visual acuity with the inventive lenses. 

Fig. 7 and 8 show modulation transfer functions of an individual with an individually 

designed lens. 

Fig. 9 shows the best corrected visual acuity with individually designed lenses according 
to the invention. 

Fig. 10 shows the age distribution of 71 patients used in a study described below in the 
example part. 

Fig. 1 1 shows a height map given by an Orbscan® true height data file. 
Fig. 12 shows average corneal wavefiront aberration coefficients. 

Fig. 13 shows a scatter plot of the spherical aberration of 71 subjects for a 6 mm diameter 
aperture. 

Fig. 14 shows a scatter plot of the spherical aberration of 71 subjects for a 4 mm diameter 
aperture. 

Fig. 15 shows a scatter plot of the spherical aberration of 71 subjects for a 5 mm diameter 
aperture. 

Fig. 16 shows a scatter plot of the spherical aberration of 71 subjects for a 7 mm diameter 
aperture. 
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Example 1 

A sample set of 10 corneal surfaces from individuals were described using 
Zernike polynomials. The sag data of the corneas was determined using the real height 
data measured with a Humphrey Atlas corneal topographer. The corneal topographer 
5 measures the height (zj) at a discrete number of points. The corneal surface can then be 
expressed as a linear combination of the first 1 5 Zernike polynomials (as described in 
Equation 1, above), where Z\ is the ith Zernike polynomial and aj is the weighting 
coefficient for this polynomial. The Zernike polynomials are a set of complete orthogonal 
polynomials defined on a unit circle. These polynomials as listed in Table 1 above and 
10 the weighting coefficients (aj) are calculated from the height data using a Grahm-Schmidt 
orthogonalization procedure. The Zernike coefficients (a*) for the 10 sample corneas are 
listed in Table 2 in mm. 
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Ann 


AOA 






rUA 


rUM 


FCZ 


_ A 

a i 


~7 ~7 4 OOQ7 


-o.oyoo4o 


-/ .ZZZOOO 


7 1^0007 


-/ .UU1 OOO 


7 oncoii 
-/ .o22b24 


-7.03713 


a2 


0.000271 


-u.uuuyoo 


U.UUUooo 


-u.uuuoiy 


U.UUU42D 


-0.000094 


-0.000236 


a3 


0.000478 


-0.000002 


-0.000847 


O.OOOyyb 


-0.000393 


0.000045 


0.000454 


a4 


0.073309 


0.083878 


0.077961 


0.078146 


0.08011 1 


0.077789 


0.079987 


a5 


-0.000316 


-0.000753 


0.000119 


0.000347 


-0.001 197 


0.00022 


-0.000071 


a6 


0.001661 


0.000411 


-0.000148 


-0.000386 


0.000314 


0.000669 


0.00079 


a/ 


0.000193 


0. 0000b 


-0.000.130 


U.UUUoz4 


-0.000161 


-0.000058 


0.000148 


a8 


0.000098 


-0.000437 


0.000146 


-0.00018 


0.000147 


0.000039 


-0.000076 


a9 


-0.000091 


-0.000168 


-0.000107 


0.000047 


-0.000181 


-0.000154 


-0.000043 


a10 


-0.000055 


0.000139 


-0.000132 


-0.000149 


0.000234 


-0.000228 


0.000244 


a11 


0.000277 


0.000394 


0.000203 


0.000305 


0.000285 


0.000315 


0.000213 


a12 


-0.000019 


-0.000105 


0.000025 


0.00007 


-0.000058 


-0.000033 


0.00009 


a13 


0.000048 


0.000032 


0.000085 


0.000017 


0.000039 


0.000059 


0.000022 


a14 


-0.000067 


0.000041 


-0.000081 


-0.000049 


0.000118 


-0.000108 


0.000127 


a15 


-0.000048 


-0.000075 


-0.000073 


-0.000019 


-0.000036 


-0.000119 


-0.000021 





FGP 


JAE 


JBH 


a1 


-7.84427 


-7.582005 


-6.890056 


a2 


-0.00056 


-0.000344 


-0.000155 


a3 


0.000347 


0.000246 


-0.000558 


a4 


0.072595 


0.075803 


0.081415 


a5 


0.000686 


-0.000388 


-0.000269 


a6 


-0.00048 


0.001688 


0.001492 


a7 


0.00014 


0.000104 


-0.000227 


a8 


-0.00025 


-0.000173 


-0.000116 


a9 


0.000092 


-0.000023 


-0.000109 


a10 


-8.2E-05 


-0.000004 


0.000065 


a11 


0.000308 


0.000309 


0.0004 


a12 


-2E-06 


-0.000115 


-0.00011 


a13 


0.000101 


-0.000042 


-0.000052 


a14 


-1.9E-05 


-0.000068 


0.00001 


a15 


0.000022 


-0.000013 


-0.000048 



Table 2: The Zernike coefficients for the 10 individual corneal surfaces in mm. 

5 These wavefront aberration coefficients can be calculated using optical design 

software such as OSLO (Sinclair Optics). Table 3 shows the results of calculating the 
wavefront aberration for subject FCM. (N.B. The normalization factor for the 
polynomials used in OSLO is different from those shown in Table 3. This difference has 
been incorporated into the coefficient values.) 
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Aberration Coefficients for FCM (OSLO) 


AO 


-0.000123 


A1 


4.5960e-07 


A2 


2.0869e-07 


A3 


-5.355e-06 


A4 


0.000551 


A5 


0.000182 


A6 


3.7296e-05 


A7 


-5.5286e-05 


A8 


0.000116 


A9 


-0.000217 


A10 


-0.000147 


A11 


-3.8151e-05 


A12 


6.1808e-05 


A13 


-3.3056e-07 


A14 


4.888e-07 


A15 


-1.8642e-06 


A16 


-0.000115 


A17 


-0.000127 



Table 3: The corneal aberration coefficients in mm calculated for subject FCM using 

OSLO( N.B. OSLO numbering order) 



Example 2 

An averaged design embodiment of the inventive lenses has been calculated 
using the average "old" cornea information provided by Pablo Artal, Murcia, Spain. This 
5 data was taken from a population sample of 16 old corneas in which all of the subjects 
had a visual acuity of 20/30 or better. The corneal surfaces were described using Zernike 
polynomials for an aperture of 2.0 mm radius (r 0 ). The polynomial coefficients were then 
used to determine the radius and asphericity values using Equations 2 and 3. 
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r 2 

R= r=-^ 7= (2) 

2(2V3a 4 -&j5a n ) 



K 2 =$£-6j5a u (3) 
r. 



Note that the asphericity constant, K, describes the surface's variation from a 
sphere (K 2 =l-e 2 ). (i.e. For a sphere K = 1 and for a parabola K = 0). (cc = K 2 -l, wherein 
cc is the conical constant) 
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Because the cornea surface has only been described for a central diameter of 4 
mm, the calculated R and K are also only accurate over the central 4 mm. A pupil size of 
4.0 mm is therefore selected for design purposes. This pupil size is reasonable for 
intraocular lens design purposes. 
5 A 22D CeeOn® 911 lens from Pharmacia Corp was chosen as a starting point for 

the averaged lens design. For the purpose of comparison, the averaged lenses were also 
designed to be 22D. (Note that other dioptres would give similar simulation results, 
provided that the spherical surfaces of the lenses is the same.) The surface information for 
the starting point eye model is summarized in Table 4. In the conical and aspherical data 
10 demonstrated in Table 4, the average conic constant CC is determined for the 10 
individual corneas of Example 1. 



Surface 
# 


Radius 
(mm) 


Thickness 
(mm) 


Aperture 
Radius (mm) 


Conic 
Constant 


Refractive 
Index 


Object 




OO 


2.272611 




1.0 


1 

(cornea) 


7.573 


3.6 


2.272611 


-0.0784 * 


1.3375 


2 (pupil) 






2.0 




1.3375 


3 




0.9 


2.0 




1.3375 


4 (lens 

i) 


11.043 


1.14 


3.0 




1.4577 


5 (lens 
2) 


-11.043 


17.2097 


3.0 




1.336 



15 Table 4: Surface data for the starting point of the averaged ("Zll") design 

* This conic constant for the "average" cornea is taken from the published works of 
Guirao and Artal 

20 

The wavefront aberration coefficients in mm for the average cornea are shown in 
column 1 of Table 5, while the coefficients in mm of the combination of the average 
cornea and the 91 1 lens are shown in column 2 of Table 5. Note that the Zll coefficient 
(al 1) of the average old cornea alone is 0.000220 mm, while the Zl 1 of this eye 
25 implanted with a 91 1 would be 0.000345 mm. 
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Average Cornea 


Average Cornea +91 1 


a1 


0.000432 


0.000670 


a2 


0.0 


0.0 


a3 


0.0 


0.0 


a4 


0.000650 


0.00101 


a5 


0.0 


0.0 


a6 


0.0 


0.0 


a7 


0.0 


0.0 


a8 


0.0 


0.0 


a9 


0.0 


0.0 


a10 


0.0 


0.0 


a11 


0.000220 


0.000345 


a12 


0.0 


0.0 


a13 


0.0 


0.0 


a14 


0.0 


0.0 


a15 


0.0 


0.0 



Table 5: Zernike coefficients in mm of the average cornea and the starting point for 
design (Average cornea + 911) 



15 The averaged lens was optimized to minimize spherical aberration, while 

maintaining a 22D focal power. The lens material remained the same as in the 22D 911 
lens (HRI silicone, the refractive index of which is 1 .45 in77 at 37°C . The resulting 
design for an equiconvex lens is The averaged lens was optimized to minimize spherical 
aberration, while maintaining a 22D focal power. The lens material remained the same as 

20 in the 22D 911 lens (HRI silicone, the refractive The resulting design for an equiconvex 
lens is provided in Table 6. The total-eye Zl 1 coefficient of the average cornea combined 
with this lens is — 2.42x 10" 7 mm(versus 0.000345 mm for the cornea plus 91 1 lens). 
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Surface # 


Radius 


Thickness 


Aperture 


Conic 


4 m Order 




(mm) 


(mm) 


Radius (mm) 


Constant 


Aspheric 
Constant 


Object 




00 


2.272611 






1 (cornea) 


7.573 


3.6 


2.272611 


-0.0784 




2 (pupil) 






2.0 






3 




0.9 


2.0 






4 (lens 1) 


10.0 


1.02 


3.0 


-2.809 


-0.000762 


5 (lens 2) 


-12.0 


17.2097 


3.0 







Surface # 


6 m Order 


Refractive 




Aspheric 


Index 




Constant 




Object 




1.0 


1 (cornea) 




1.3375 


2 (pupil) 




1.3375 


3 




1.3375 


4 (lens 1) 


-1.805e-05 


1.4577 


5 (lens 2) 




1.336 



5 

Table 6: Surface data for the starting point of the averaged lens design 

The corneas of the 10 test subjects were combined in an optical system with both 
the 911 and the averaged lenses. The resulting total-eye Zl 1 coefficients are shown in 

10 Fig. 1 . As demonstrated, in Fig. 1, in each case, the absolute value of the Zl 1 coefficient 
was less when the Zl 1 lens was implanted. Because subjects CGR and FCZ have 
relatively low levels of corneal spherical aberration to begin with, the total-eye spherical 
aberration is overcorrected in these two cases. As a result, the sign of the total spherical 
aberration is noticeably reversed in these two cases, and the amount of spherical 

15 aberration is still significant. In every other case, the spherical aberration of the total eye 
would be essentially 0 after the implantation of a Zl 1 lens. The visual acuity of each of 
the 10 test subjects were calculated according to standard methods described in "Visual 
acuity modeling using optical raytracing of schematic eyes", Greivenkamp et al., 
American journal of ophthalmology, 120(2), 227-240, (1995), for the implantation of 

20 both a 22D 91 1 lens and a 22D averaged "Zl 1" lens. The square wave responses were 

calculated using OSLO™ and a software module was written in Matlab™ to calculate the 
visual acuity following the above method. The resulting visual acuities are shown in Fig. 
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2. Out of the 10 cases investigated and shown in Fig. 2, eight subjects had better vision 
when implanted with the averaged lens according to the present invention. In the cases 
where the visual acuity decreased their Snellen distance increased by less than 1 ft which 
would not show up in visual acuity testing. 

To be able to assess the optical quality difference between a CeeOn® 91 1 A and 
averaged lenses according to the present invention, a physical model of the average 
cornea was designed and manufactured. It is a convex-piano lens of PMMA with an 
aspheric front surface having a value of 0.000218 for Zernike coefficient al 1. This value 
is almost equal to the value of the calculated average cornea: 0.000220. With the PMMA 
model cornea MTF measurements were performed on an optical bench in a model eye 
with the "averaged" Zl 1 lenses and CeeOn® 91 1 A lenses. Modulation Transfer Function 
(MTF) measurements are a widely accepted method of quantifying image quality. 
Through focus MTF measurements at 50c/mm and a frequency MTF curves focussed at 
50c/mm, in both cases with a 3mm pupil are shown in the Fig. 3 and Fig. 4, respectively 
for lenses with an optical power of 20D. The width of the through focus MTF at 0.2 MTF 
units is a measure for the depth of focus and is equal for both lenses. The MTF curve 
focussed at 50c/mm for "averaged" ZI 1 lenses is almost diffraction limited and is better 
than that for CeeOn 91 1 A lenses. . 

The astigmatism of the cornea and the defocus of the system can be corrected by 
adjusting the Zernike coefficients of the cornea model and the focal position of the 
system. When this is done and the procedure for calculating visual acuity is repeated the 
results in Fig. 6 are obtained. They represent a modeled best corrected visual acuity. We 
now see that, in all cases, after correction for astigmatism and defocus (as in reality 
would be done with spectacles) the inventive averaged lens produces a higher best 
corrected visual acuity than the 91 1 lens of the same dioptre. 
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Example 3 

Individually Designed Lenses: 

As a potential further improvement upon the averaged lens ("Zl 1 lenses"), an 
individualized lens ("II 1 lenses 55 ) was designed for each of four subject corneas 
employing the same design principles as demonstrated in Example 2. The individual 
lenses were designed such that the Zl 1 of the lens balances the Zl 1 of the individual 
cornea. The total-eye Zl 1 coefficients for the II 1 lenses are shown in Table 7, together 



Subject 


911 


averaged 


individual 


CGR 


0.000107 


-0.000101 


-0.000004 


FCZ 


0.000105 


-0.000100 


-0.000005 


JAE 


0.000195 


-0.000016 


-0.000012 


JBH 


0.000238 


0.000037 


-0.000019 



Table 7: The Zl 1 coefficients in mm of the model eyes with the 91 1, Zl 1 and 111 lenses 



with the corresponding coefficients for the 911 and the averaged lenses. Furthermore, for 
each of the 91 1, Zl 1 (averaged), and 111 (individual) lenses, the MTF curve at best 
focussed at 50 c/mm and the through focus MTF at 50c/mm for subject JAE are plotted 
below in Fig. 7 and 8. From Fig. 7 and 8, it is seen that the MTF at 50 c/mm of eyes 
implanted with the Zl 1 and II 1 lenses is higher than the MTF of the same eyes implanted 
with 91 1 lenses. It can also be seen that the through focus MTF of all of the lenses is 
satisfactory. The Zl 1 has as much depth of focus as the 911. However, it is also 
interesting to note that the 111 does not provide a significant improvement in either MTF 
or through focus MTF, relative to the Zl 1 lens. 

The visual acuities of the subjects with individualized lenses have also been calculated. 
Fig. 9 compares these acuities with the visual acuity calculated for the 911 and Zl 1 
lenses. 

From Fig. 9, we see that for all 4 subjects, visual acuity is better for both the Zl 1 and 111 
lenses than it is for the 91 1 lens. We also see that the results with the Zl 1 and II 1 lenses 
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do not differ significantly - the average cornea is relatively accurate for each of the 4 test 
subjects. 

Example 4 

5 

The design of an ophthalmic lens, which is adapted to reduce the spherical 
aberration of an average cornea obtained from a group of people will be described in 
detail here below. The lens will be called the Zl 1 lens because it compensates for the 
normalized 1 1 th Zernike term describing spherical aberration of the corneas. It was 
10 decided to use a population of potential recipients of the Zl 1 lens, namely cataract 
patients. 

Description of the Population: 

15 The population included 71 cataract patients from St. Erik's eye hospital in 

Stockholm, Sweden. These patients were of ages ranging from 35 to 94 years (as of April 
12, 2000). The average age of our population was 73.67 years. A histogram of the age of 
the population is shown in Figure 10. 

The corneas of the 71 subjects were measured using an Orbscan® (Orbtek, Salt 

20 Lake City) device. Orbscan® is a scanning slit -based, corneal and anterior segment 

topographer that measures both surfaces of the cornea, as well as the anterior lens surface 
and the iris. Each surface can be displayed as maps of elevation, inclination, curvature, 
and power. 

25 Fitting Algorithm: 

The corneal elevation height data (the Cartesian locations of points on the surface 
of the cornea) for the anterior surface was obtained using the Orbscan®, and used as raw 
data for the determination of the optical properties of the cornea. The height data from an 
30 example Orbscan® file is represented in Figure 1 1 . 
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The Cartesian co-ordinates representing the elevation height data are transformed 
to polar co-ordinates (x,y,z — ► r,0,z). In order to describe the surface, this data is then fit 
to a series of polynomials as described in Equation lb. The coefficients (a's), or 
weighting factors, for each polynomial are determined by the fitting procedure resulting 
in a complete description of the surface. The polynomials (Zj) used are the normalized 
Zernike polynomials. 



10 These polynomials are special because they are orthonormal over a continuous 

unit circle. They are commonly used to describe wavefront aberrations in the field of 
optics. Corneal topographers measure the elevation heights at a discrete number of points. 
The Zernike polynomials are not orthogonal over a discrete set of points. However, 
applying an orthogonalization procedure, termed Gram-Schmidt orthogonalization, to the 

15 height data, allows the data to be fit in terms of Zernike polynomials maintaining the 

advantages of an orthogonal fit. Sixty-six coefficients (a's) were used to fit the height data 
provided by the Orbscan® software. A Matlab™ algorithm was used in the fitting 
procedure. The radius and asphericity value can be approximated from the Zemike 
coefficients (Equations 2b and 3b) and the conic constant of the surface is simply K 2 - 1 

20 (from this we know that for a sphere K 2 = 1). The fitting procedure is well described in a 
number of references. Four different articles are refereed to here: "Wavefront fitting with 
discrete orthogonal polynomials in a unit radius circle", Daniel Malacara, Optical 
Engineering, June 1990, Vol. 29 No. 6, "Representation of videokeratoscopic height data 
with Zemike polynomials", J. Schwiegerling, J.Greivenkamp and J.Miller, JOSA A, 

25 October 1995, Vol. 12 No. 10, "Wavefront interpretation with Zemike polynomials" 

J.W.Wang and D.E. Silva, Applied Optics, May 1980, Vol. 19, No. 9 and "Corneal wave 
aberration from videokeratography: accuracy and limitations of the procedure", Antonio 
Guirao and Pablo Artal, J Opt Soc Am A Opt Image Sci Vis Jun 2000, Vol. 17(6):955-65. 



R = r Z*" (2b ) 



2(2j3a 4 -&j5a n ) 
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K m =%6^5a u (3b) 

Calculation of Wavefront Aberration: 

Knowing the shape of the anterior corneal surface (Zernike coefficients described 
above as a's) 5 it is possible to determine the wavefront aberration contributed by this 
surface using a raytrace procedure. This is described in for example "Corneal wave 
aberration from videokeratography: accuracy and limitations of the procedure", Antonio 
Guirao and Pablo Artal, J Opt Soc Am A Opt Image Sci Vis Jun 2000, Vol. 17(6):955-65. 

Results: 

Average Corneal Spherical Aberration and Shape: 

The 71 subjects were evaluated using the criteria described above for a 6 mm 
aperture. The wavefront aberration of each subject was determined after fitting the 
surface elevation with Zernike polynomials. Figure 12 shows that average and standard 
deviation of each Zernike term (normalized format). The error bars represent ± 1 standard 
deviation. There are three aberrations that are significantly different from zero on average 
in our population. These are astigmatism (A5), coma (A9) and spherical aberration (All). 
Spherical aberration is the only rotationaliy symmetric aberration, making it the only 
aberration that can be corrected with a rotationaliy symmetric IOL. 

Figure 1 3 shows a scatter plot of the value of the Zernike coefficient (OSLO 
format) representing spherical aberration for each of the 71 subjects before cataract 
surgery. The solid line in the middle represents the average spherical aberration, while the 
dotted lines represent +1 and —1 standard deviation. Table 8 lists the average, standard 
deviation, maximum and minimum values for the radius, aspheric constant, spherical 
aberration and root mean square error. 
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Average Value 


Standard Deviation 


Maximum 


Minimum 


R(mm) 


7.575 


0.333 


8.710 


7.072 


Ksq 


0.927 


0.407 


2.563 


0.0152 


SA coefficient 
OSLO format (in 
mm) 


0.000604 


0.000448 


0.002003 


-0.000616 


RMSE 


0.000055 


0.00000482 


0.000069 


0.000045 



Table 8: the average, standard deviation, maximum and minimum values for the radius, 
aspheric constant, spherical aberration and root mean square error for a 6mm aperture. 

Tables 9,10 and 1 1 below show the corresponding results for aperture sizes of 4,5 
and 7 mm respectively. Figure 14,15 and 16 are the corresponding scatter plots. 





Average Value 


Standard Deviation 


Maximum 


Minimum 


R 


7.56292 


0.320526 


8.688542 


7.067694 


Ksq 


0.988208 


0.437429 


2.33501 


-0.051091 


SA(A11 
in mm) 


0.000139 


0.000103 


0.00041 


-0.000141 


RMSE 


4.52E-05 


4E-06 


0.000054 


0.000036 
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Table 9: The average, standard deviation, maximum and minimum values for the radius, 
aspheric constant, spherical aberration and root mean square error using an aperture 
diameter of 4 mm. 
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Average Value 


Standard Deviation 


Maximum 


Minimum 


R 


7.55263 


0.320447 


8.714704 


7.09099 


Ksq 


0.945693 


0364066 


2.045412 


0.044609 


SA 

(All in mm) 


0.00031189 


0.000208 


0.000793 


-0.000276 


RMSE 


4.7E-05 


4.02E-06 


0.000057 


0.000037 



Table 10: The average, standard deviation, maximum and minimum values for the radius, 
aspheric constant, spherical aberration and root mean square error using an aperture 
diameter of 5 mm. 



5 





Average Value 


Standard Deviation 


Maximum 


Minimum 


R 


7.550226 


0.336632 


8.679712 


7.040997 


Ksq 


0.898344 


0.416806 


2.655164 


-0.04731 


SA 

(All in mm) 


0.001058 


0.000864 


0.003847 


-0.001319 


RMSE 


7.58E-05 


L02E-05 


0.000112 


0.000057 



Table 1 1 : The average, standard deviation, maximum and minimum values for the radius, 
aspheric constant, spherical aberration and root mean square error using an aperture 
diameter of 7 mm. 
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Design Cornea: 

One model cornea was designed and each Zl 1 lens power was designed using this 
cornea. The cornea was designed so that it had a spherical aberration that is the same as 
15 the average calculated for the population. The design cornea radii and aspheric constants 
are listed in Table 12 for different aperture sizes. In every case, the radius of curvature 
was taken to be the average radius determined from the Zernike fit data. The aspheric 
constant was varied until the spherical aberration value of the model cornea was equal to 
the average spherical aberration value for the population. 

20 
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Mpenure size ^rrirn; 


PvcaUIUS ^1 III II y 


iconic oonsiani 
(OSLO value, K 2 -1) 


i i ^oemciem (rnrnj 


4 


7.563 


-0.0505 


0.000139 


5 


7.553 


-0.1034 


0.000312 


6 


7.575 


-0.14135 


0.000604 


7 


7.55 


-0.1810 


0.001058 



Table 12: The design cornea radii and aspheric constants for aperture diameters of 4, 5, 6, 
and 7 mm. 

As discussed previously, the 6 mm aperture diameter values are used for the 
design cornea. This choice enables us to design the Zl 1 lens so that it has no spherical 
aberration (when measured in a system with this cornea) over a 5.1 mm lens diameter. 
The OSLO surface listing for the Zl 1 design cornea is listed in Table 13. The refractive 
index of the cornea is the keratometry index of 1 .3375. 

These values were calculated for a model cornea having a single surface with a 
refractive index of the cornea of 1.3375. It is possible to use optically equivalent model 
formats of the cornea without departing from the scope of the invention. For example 
multiple surface corneas or corneas with different refractive indices could be used. 



Surface # 


Radius 


Thickness 


Aperture 


Conic 


Refractive 




(mm) 


(mm) 


Radius 


Constant 


index 








(mm) 


(cc) 




Object 




1.0000e+20 


1.0000e+14 




1.0 


1 (cornea) 


7.575000 


3.600000 


3.000003 


-0.141350 


1.3375 


2 (pupil) 






2. 640233 




1.3375 


3 




0. 900000 


2 . 64023 




1. 3375 


4 




25. 519444 


2.550292 




1. 3375 


5 






2.2444e-05 




1.3375 



Table 13: OSLO surface listing for the Zll design cornea. 
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Lens Design: 

Each Zl 1 lens was designed to balance the spherical aberration of the design 
cornea. The starting point for design was the CeeOn Edge© 91 1 lens described in U.S. 
5 5^444^106 of the same power, with modified edge and center thickness. The lens was then 
placed 4.5 mm from the anterior corneal surface. The distance from the anterior corneal 
surface is not that critical and could be varied within reasonable limits. The surface 
information for the starting point eye model for the 22 D lens design process is listed in 
Table 14. The anterior surface of the lens was described using the formula shown in 

10 Equation 4. The variables cc, ad and ae were modified to minimize the spherical 

aberration. The variables are determined for an aperture size of 5.1 mm and the surface is 
extrapolated from these values to the optical aperture size of 6 mm. The resulting 22D 
Zl 1 eye model is listed in Table 15. The anterior surface of this 22D lens has been 
modified in such a way that the spherical aberration of the system (cornea+lens) is now 

15 approximately equal to 0. The wavefront aberration coefficients as calculated by OSLO 
for the CeeOn Edge 911 22D lens eye model and the 22D Zl 1 lens eye model are listed 
below in Table 16. Note that the coefficient representing spherical aberration for the 
starting point eye model is 0.001005 mm for a 6mm diameter aperture placed at the 
cornea, while the same coefficient for the eye model with the designed Zl 1 lens is - 

20 1 .3399e-06 mm. The same process as described above for a 22D lens can similarly be 
performed for any other lens power. 



z = 




(4) 
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Surface # 


Radius 


Thickness 


Aperture 


Conic 


Refractive 




(mm) 


(mm) 


Radius 


Constant 


index 








(mm) 


(cc) 




Ob j ect 




1 . uuuue + ^i) 


1 . UUUUe-rl4 




1 . 0 


1 (cornea) 


7.575 


3.600000 


3.000003 


-0.14135 


1.3375 


2 (pupil) 






2. 640233 




1.336 


3 




0.900000 


2 . 64023 




1.336 


4 (lens) 


11.043 


1 .164 


2.550191 




1 .458 


5 (lens) 


-11.043 


17 . 1512 


2.420989 




1.336 


6 (image) 


0.0 


-0.417847 


0. 058997 







Table 14: Surface data for the starting point averaged eye model and a 22D lens 
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Cn ■>- fa 
# 


(mm) 


TYi icknes s 
(mm) 


Aperture 
Radius 
(mm) 


Conic 
Constant 
(cc) 


4 th Order 
aspheric 
Constant 


Object 




1.0e+20 


1.00e+14 






1 

(cornea) 


7.575 


3. 60 


3. 00 


-0. 14135 




2 

(pupil) 






2 , 64 






3 




0. 90 


2. 64 






4 (lens ) 


11.043 


1. 164 


2.55 


-1.03613 


-0. 000944 


5 (lens) 


11.043 


17.1512 


2.42 






6 (image ) 






1.59e-05 







Surface 
t 


6th Order 
aspheric 
constant 


Ref ractiv 
e index 


Object 




1.0 


1 

(cornea) 




1. 3375 


2 

(pupil) 




1.336 


3 




1.336 


4 ( lens ) 


-1.37e-05 


1 . 458 


5 (lens) 




1 .336 


6 (image) 







Table 15: Surface data for the averaged eye model and the final 22D Zl 1 lens 
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uoemcient 


Average corneal- 22D 911 


Average cornea + llu Zl 1 


al 


A AAAA^C") 
-0. 000962 


1 on/:« a/t 

-1.89oe-06 


_ 

az 


A A 
0.0 


A A 

0.0 


a3 


0.0 


0.0 


a4 


2.3101e-05 


-3.9504e-06 


a5 


0.0 


0.0 


a6 


0.0 


0.0 


a7 


0.0 


0.0 


a8 


0.0 


0.0 


a9 


0.00105 


-1.3399e-06 


alO 


0.0 


0.0 


all 


0.0 


0.0 


al2 


0.0 


0.0 


al3 


0.0 


0.0 


al4 


0.0 


0.0 


al5 


0.0 


0.0 



Table 16: Zernike coefficients (OSLO format) for the average cornea and a 22D 911 lens 
and the average cornea and the 22D Zl 1 lens 

5 

The optical form chosen for the new Zl 1 design is an equiconvex lens made from 
a silicone with refractive index of 1.458. The spherical aberration of an average cornea is 
balanced by the Zl 1 lens yielding a system without spherical aberration. The front 
surface of the lens is modified such that the optical path lengths of all on-axis rays within 

10 the design aperture are the same producing a point focus. This feature can be achieved 
with many lens forms. The Zl 1 lens could therefore be designed on a convex-piano, 
plano-convex, non-equiconvex lens or any other design yielding a positive lens. The Zl 1 
concept could also be extended in order to encompass a negative lens used to correct the 
refractive errors of the eye. The front surface or back surface could also be modified to 

15 produce the needed change in optical path difference that neutralizes the spherical 
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aberration. There are therefore many possible designs that would achieve the goals of the 
Zll lens design. 

5 



